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ABSTRACT 

Gyrochronology and recent theoretical findings are used to reduce the number of input 
parameters of differential rotation models. This eventually leads, after having fixed our 
turbulence model parameters, to a theoretical prediction for the surface differential 
rotation as a function of only two stellar parameters - surface temperature and ro- 
tation period - that can be defined observationally. An analytical approximation for 
this function is suggested. The tendency for the differential rotation to increase with 
temperature is confirmed. The increase is much steeper for late F-stars compared to 
G- and K-dwarfs. Slow and fast rotation regimes for internal stellar rotation are iden- 
tified. A star attains its maximum differential rotation at rotation rates intermediate 
between these two regimes. The amplitude of the meridional flow increases with sur- 
face temperature and rotation rate. The structure of the flow changes considerably 
between cases of slow and fast rotation. The flow in rapid rotators is concentrated in 
the boundary layers near the top and bottom of the convection zone with very weak 
circulation in between. 
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1 INTRODUCTION 

Rotation of the sun and other stars with external convec- 
tion zones is not uniform. The rotation period depends on 
latitud e and typically incr e ases from the equator to th e 
poles l|Donahu'e et aLMl996l : ICollier Cameron! |2003l . \200H ). 
The phenomenon of differential rotation is closely related 
to magnetic activity and has long been a subject of theo- 
retical stud ies. The currently leading theoretical concept pi- 
oneered by iLebedinskiil (Il94ll ) explains differential rotation 
by the interaction between convection and rotation. Con- 
vective motions in a rotating star are disturbed by Coriolis 
force. Back reaction disturbs rotation to make it not uni- 
form. 

Theoretical models rep roduce the observed rotation 
of the sun quite closely dKitchatinov fc Riidigerl 120051 : 
iMiesch et alll2006t ). In the case of the sun, however, the 
differential rotation theory has to explain what is already 
known from helioseismology. The detailed h elioseismolog- 
ical p icture of the inte rnal solar rotation (| Wilson et al.l 
1 19971 : ISchou et all Il998l ) left no opportunity for theoret- 
ical predictions. The theory then has to test its predic- 
tive ability with other stars, and this is very tempting 
to do in view of the rapid development of asteroseismol- 
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ogy (|Christensen-Dalsgaard fc Hou dck 201(J). Some com- 
putations of differential rotation as a function of stel- 
lar mass and rotation rate have been already attempted 
(IKitchatinov fc Rudigerlll99il and were even to some extent 
supported by observations I Barnes et al. 1 120051 ). A system- 
atic study of stellar differential rotation was, however, not 
possible because the stellar parameters on which it depends 
are too many. The resulting differential rotation depends on 
mass, chemical composition (metallicity), rotation period, 
and age of a star. Exploring four-dimensional parametric 
space is an unbearable task. 

The situation has changed recently. First, the develop- 
ment of a new branch of astronomy named Gyrochronology 
has led to the establishment of a functio nal relation between 
stellar mass, age, and rate of rotation llBarnesll200l 120071 . 
l2010l : ICollier Cameron et al.l 12009] ; iMeibom et al.1 120091 ). 
Second, theoretical modeling shows that main-sequence 
stars of different mass and metallicity have (almost) the 
same diff erential rotation when their su rface temperatures 
coincide (jKitchatinov fc Olemskovl liuTll . hereafter KOll). 
In other words, the dependencies on mass and metallicity 
can be combined into a common dependence on effective 
temperature. These two developments reduce the number of 
independent parameters on which differential rotation de- 
pends by two. The dependence on the remaining two pa- 
rameters can well be explored numerically. 
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In our preceding publication (KOll), a new mean- 
field model for differential rotation was presented. Perfor- 
mance of the model was tested by applying it to the sun 
and individual stars including two rapid rotators (AB Dor 
and LQ Hya) and two moderate rotators (e Eri and 
Ceti) whose rotation has been studied well observation- 



ally dDonati fc Collier Cameron! Il997l ; iKovari et ail 12004 
ICroll et al.ll2006l ; IWalker et all 120071 ). The temperature de- 
pendence of differen tial rotation of young stars observed by 
iBarnes et ail (2005) was modeled. This new paper reports 
and summarizes the results of extensive computations of 
stellar differential rotation. We find that after a solar-like 
star arrives at the main sequence, its differential rotation 
initially increases and then decreases as the star ages. The 
rotation period, at which the differential rotation is maxi- 
mum, decreases with stellar mass. The dependence on ro- 
tation rate is, however, mild. The differential rotation de- 
pends more strongly on surface temperature, the hotter the 
star, the larger it is. In this paper, we suggest an analyti- 
cal approximation for the numerically defined dependence of 
differential rotation on rotation period and surface temper- 
ature. The theoretical prediction can be tested by observa- 
tions if only rotation period, surface temperature and the 
surface differential rotation are simultaneously measured. 
Dependence of the meridional flow amplitude and struc- 
ture on rotation rate and temperature is also studied. The 
next Section [2] describes our model and method. Section 
presents and discusses the results, which are summarized in 
Section 0] 



2 MODEL AND METHOD 
2.1 The model 

Our model is based on the mean-field theory of differential 
rotation. This means that the velocity field (v) is split in its 
mean (V) and fluctuating (it) parts, 



v = V + u, V={v), (u)=0, 



(1) 



and equations for the mean flow, V, are formulated. The 
equations are then solved numerically. The angular brack- 
ets in (fTJ) signify an averaging procedure that smoothes out 
random fluctuations of turbulent fields. The steady equation 
for the mean flow, 



(V- 



V) V + -VP 

p 



1 



(2) 



involves the effects of turbulence via the Reynolds stress 
tensor, 



Rij — —p( u i u j)- 



(3) 



In equation p, P, and g are the mean density, pressure, 
and gravity respectively. Reynolds stress (O includes two 
parts representing the turbulent viscosity and non-viscous 
flux of momentum, 



Rij — ~PQii + pNijkl—. 

ori 



(4) 



where the summation convention over the repeated sub- 
scripts is assumed, Mijki is the eddy viscosity tensor, and 
Qij represents the non-viscous part of the stress (cf. Ap- 
pendix of KOll for detailed expressions for both parts of 



the Reynolds stress @). The presence of non- viscous flux of 
angular mome ntum in rotati ng turbulent fluids was called 
the 'A-effect ' i|Riidigerl 1989). This effect is the principal 
driver of differential rotation in our model. 

The direction of the angular momentum transport by 
the A-effect is controlled by two dimensionless parameters, 
V {0) and H w , of the Q A -tensor 



Q'i 



v T Q k fi (V (0) (fT) (fi£j k i + fjSiki) 



+ 



(5) 



where SI is the angular velocity vector, r is the radial unit 
vector, 



r£ 2 g dS 
15c p dr 



(6) 



is the eddy viscosity parameter, S is the specific entropy, g 
is gravity, c p is the specific heat at constant pressure, £ — 
a MLT f/p is the mixing length, H p is the pressure scale height, 

SI* = 2tQ (7) 

is the Coriolis number, 

1/3 
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(8) 



is the convective turnover time, 6F = L(4irr 2 ) 
the 'residual' heat flux that convection has to transport, and 
P rad is the radiative heat flux. 

Coriolis number measures the intensity of inter- 
action between convection and rotation. Its value defines 
whether the turbulent eddies are living long enough for ro- 
tation to influence them considerably. The Coriolis num- 
ber, ST = 47TRO" 1 , is reciprocal to the Rossby number 
Ro = ProtT -1 ; Prot is the rotation period. Explicit ex- 
pressions for th e functions V (0) (Q*) a nd H (1) (Q.*) of (5J 
were derived bv iKitchatinov fc Riidieerl (1995,2005). In the 
slow rotation limit, fl* < 1, H^(Q*) is small compared 
to V (0) (Q*) and V (0) (fT) is negative. This means that the 
non-diffusive flux of angular momentum is close to the in- 
ward radial direction. In the opposite limit of rapid rotation, 
fT > 1, y (0) (fT) is small compared to H {1) (SI*). Angular 
momentum flux is parallel to the rotation axis and points 
to the equatorial plane in this case. As the Coriolis number 
increases from small to large values, the angular momentum 
transport by the A-effect changes smoothly from radial in- 
ward to equatorward and parallel to the rotation axis. This 
picture generally agrees with 3D numerical simulation s of 
iKapyla. fc Brandenburg (|2008h and lKapyla et alj (120111 ). 

The differential rotation model of this paper is identical 
to that of KOll. A detailed formulation of the model can 
be found in that paper and in the references given there. We 
just outline the model here. 

The azimuthal component of Eq. ([2]) in spherical coor- 
dinates (r, 6, cj>) gives the mean-field equation for the angular 
momentum balance, 



div (pr sin d(u^u) + pr 2 sin 2 6 SlV m ) = 0, 



(9) 



where V" 1 is the meridional flow velocity. Upon substitution 
of an explicit expression for Reynolds stress in (O, an equa- 
tion for angular velocity can be obtained, which involves the 
meridional flow and is thus not closed. 
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The origin of meridional flow can be seen from the equa- 
tion, which can be obtained as the azimuthal component of 
the curled equation ([2]), 



V{V m ) = sinf 



an 2 

8z 



c p r 88 ' 



(10) 



where 8/8z = cos 88/dr—r^ 1 sin 88/ 89 is the spatial deriva- 
tive along the rotation axis. The left part of this equation 
describes the viscous drag on the meridional flow 



drj \pdri V 



dV„ 

8r„ 



(11) 



(see Appendix in KOll for the expression for D(V m ) in 
spherical coordinates) . 

The right part of the meridional flow equation (|10[) in- 
cludes two sources of meridional flow: the centrifugal driv- 
ing of the flow by non-conservative part of centrifugal force 
(known also as the 'gyroscopic pumping') and the baro- 
clinic driver due to the meridional gradient of entropy (ther- 
mal wind). Each term on the right side of (|10[1 is large 
compared to the left side. The two terms on the right, 
therefore, almost balance each other, this state is known 
as the Taylor-Proudman balance. The meridional flow re- 
sults from a slight deviation from the balance. The merid- 
ional flow attains its largest velocities near the boundaries of 
the convection zone because the stress-free boundary condi- 
tions are not compatible with the Taylor-Proudman balance. 
Therefore, the Taylor-Proudman balance has to be violated 
near the stress-free boundaries where the source of merid- 
ional flow in the right side of (1101) attains its largest va lues 
l|Kitchatinov fc Rudigerlll999l ; iMiesch fc Hindmanll201lf ). 
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Direction of the convective heat flux, 

nv m ^8 



(12) 



depends on the structure of the eddy conductivity tensor 



XH = xSij + Xll - 



(13) 



The rotationally induced anisotropy of the conductivity (|13[) 
results in a deviation of the heat flux (|12[) from a radial 
direction to the poles, which in turn produces the depen- 
dence of entropy on latitude. This 'differential temperature' 
is, therefore, produced not o nly by the d ependence of ther- 
mal conductivity on latitude (Weiss 196 5]) but mainl y by the 
anisotropy of the conductivity ( Riidiger et al.ll2005T ). 

The model solves the system of three joint equations 
governing distributions of angular velocity, meridional flow, 
and entropy inside the convection zone. These steady equa- 
tions are solved numerically by the relaxation methocj]. 

It should be noted that the model relies on a uni- 
form theoretical basis: all the model needs, including the 
A-effect, eddy viscosities and thermal conductivities was de- 
rived within the same quasi-linear approximation of turbu- 
lence theory. In particular, the eddy conductivities (|13[) were 
derived as fully nonlinear functions of the Coriolis number: 



X = X T <Kfi*)> X|| 



(14) 

l|Kitchatinov et all |l994). In this equation, \t ls l ne 
isotropic eddy conductivity which would take place under 

1 The numerical solver can be obtained freely on request. 



actual sources of turbulence (actual superadiabaticity) but 
in non-rotating fluid. The mixing-length relation 



r£ 2 g 8S 
' 12c p 8r ' 



(15) 



defines the conductivity in terms of the entropy gradient. 
This gradient is controlled by the solution of the model equa- 
tions. Similarly, the eddy viscosity tensor, 



+ V2 ( On — ttt: \-0ji 



+ Oik-^~ +o jk -^- +d k i——) 



^-^T^ ^ , (16) 

involves the rotationally induced anisotropy and quenching 
via dependence of the viscosity coefficients 



v n = i/ T 0„(n*) n 



(17) 



on the Coriolis number @. Explicit expressions for 
the viscosity quenching functions, <f> n (Q.*), are given in 
iKitchatinov et al. ( 1994) and the background viscosity v T 
is defined in terms of the entropy gradient in (|6]). 

In other words, all the eddy transport coefficients are 
not prescribed but computed. Such an approach helps 
to avoid arbitrariness in prescribing uncertain parameters. 
Still, two parameters remained for some time uncertain. 
First, the A-effect depends - apart from the well known 
(nearly adiabatic) density stratification - also on turbu- 
lence anisotropy, which is uncertain. The contribution of 
anisotr opy was put to zero in the former version of our 
model l|Kitchatinov fc Ri idiger 199^) but then fixed by com- 
parison with 3D f-plane simu lations of thermal convection 
IIKitchatinov fc Riidigerll2005l ) and never varied since then. 
Second, the c x parameter of (|14p is uncertain (c x = 1 
in Kitchatinov & Riidiger 1999). It was found in KOll 
that c x = 1.5 gives close agreement with helioseismology. 
Since then, c Y = 1.5 is used for stellar applications also 
ijKiiker et al.ll201lf ) and not varied any more. It is c x = 1.5 
in this paper also. 

After the two uncertain parameters were fixed, no fur- 
ther tuning is possible in the mean-field model for differ- 
ential rotation. In this sense, the model does not have free 
parameters. 



2.2 Reducing the number of input parameters 

Input parameters of the model include the rotation period 
(Prot) and a set of parameters, all of which can be taken 
from a structure model of a (non -rotating ) star. We use the 
EZ code of stellar evolution by iPaxtonl (120041 ) to specify 
the structure of a main-sequence star of given mass, age, 
and metallicity. The four basic input parameters - rotation 
period, mass, age, and metallicity - are, however, not inde- 
pendent. 

Figure [1] shows the surface differential rotation com- 
puted with our model for stars of different mass and metal- 
licity. All stars are at the age of 1 Gyr. A rotation period 
Prot = 10 days was prescribed. The dependence on metallic- 
ity for stars of given mass is quite pronounced. However, the 
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Figure 1. Surface differential rotation as function of stellar mass 
(a) and surface temperature (b) after our model. Lines of different 
styles correspond to different metallicities Z. Computations were 
performed for 1 Gyr old stars with P ro t = 10 days. 



rapid ly rotating stars of the so-called 'C-sequence' {Barnes] 
l2007h . The structure of a star, however, changes little over 
its C-sequence life. Therefore, the relation (|18p can still be 
used for our purposes even for very young stars. We used the 
tables of colour-temperature relati ons and the interpolation 
code of lVandenBerg &: Cleml {2003) to compute B — V colour 
of Eq. (fT8l). 

In this paper we fix the metallicity Z = 0.02 and com- 
pute differential rotation for main-sequence dwarfs of differ- 
ent masses and ages. The computations cover the mass range 
from 0.5Mq to I.SMq. The stellar masses were spanned by 
0.05M© in the range of 0.5 sj M/M Q < 1.1. The mass de- 
pendence of differential rotation steepens for more massive 
stars. Accordingly, the span was reduced to O.O25M0 for 
the range of 1.1 ^ M/Mq ^ 1.3. For a given mass, structure 
models were generated for a sequence of ages. Then, rotation 
periods of Eq. (|18p were estimated for every age and corre- 
sponding structure, and a sequence of differential rotation 
models was computed. This gives the differential rotation as 
a table-given function of rotation period and temperature. 
The surface differential rotation is the n approximate d by a 
neural network with one hidden layer {Conw ay 1998): 

N 

An = W + ^2\V n tanh (w„ + wZT eB + W^-Prot) , (19) 

n = l 

where AQ = Q cq — Q po ic is the angular velocity difference 
between the equator and pole. The convenience of this ap- 
proximation is related to i) its universality: any continuous 
function can be represented by one hidden layer network 
to an arbitrary accuracy given sufficiently large N, and ii) 
availability of reliable numerical routines for tuning the pa- 
rameters of the network (|19[) . Using the regression formula 
(|19|l . a theoretical prediction for the surface differential ro- 
tation of a star can be inferred from an observed rotation 
period and effective temperature. 



dependence almost disappears when differential rotation is 
considered as a function of surface temperature. The same 
scaling was found in KOll for ZAMS stars with P ro t = 1 day. 
We conclude that the values of surface differential rotation 
computed for stars of different masses and metallicities but 
fixed rotation period lie on a single line when plotted as 
a function of the effective temperature (T e g). A number of 
other computations were done to confirm this temperature- 
scaling for stars with not too low metallicities (Population I 
stars). The scaling reduces the number of input parameters 
by one: the mass and metallicity dependencies are combined 
into a common dependence on temperature. 

Further reduction is achieved by using Gyrochronology. 
The rotation rate of solar- and late-type stars is a function of 
Teg an d age (t). Gyrochronology generalizes the ISkumanicbl 
{19721 ) law, P ro t oc t ' 5 , by specifying the proportionality 
constant in this law as a function of T e ff or other equivalen t 
parameter. We use the empirical relation of iBarne si {20071 ). 

P ro t = at n (B — V — c) b days, (18) 

with n = 0.519, a = 0.773, b = 0.601, c = 0.4, to specify the 
rotation period of a star of given B — V colour and age t 
measured in Myr. The relation does not apply to very young 



3 RESULTS AND DISCUSSION 

Figures [2] and [3] show the computed differential rotation as a 
function of the rotation period for a variety of stellar masses. 
Lines of these plots were computed for a certain initial mass 
but not for a given stellar structure. The structure changes 
as a star ages and rotation period increases. 

Dwarf stars do not spin-down beyond a certain maxi- 
m um rotation perio d depending on spectral type (see Fig. 1 
of lRengaraianll9 84). The dependence can be approximated 
by the linear function of the B — V colour, 

P™ = 109 (B — V) — 43 day. (20) 

Lines of Fig. [2] and [3] are terminated shortly beyond the 
maximum rotation period of Eq. (|20[) . 

The dependence of differential rotation on rotation rate 
is rather mild. Figures [2] and [3] show much stronger depen- 
dence on stellar mass, with more massive and hotter stars 
having larger differential rotation. The temperature depen- 
dence steepens as temperature increases. The dependence is 
in at least qualitative agreement with differential r o tation 
measurements by Doppler imaging of 1 arncs et all {20051 ) 
(for quantitative comparison of obs erved and co mputed ro- 
tation of young stars see Fig. 4 of iKitchatino-Jl 20111. The 
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Figure 2. Surface differential rotation as a function of the rota- 
tion period for stars of subsolar masses. 
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Figure 3. Surface differential rotation as a function of the rota- 
tion period for stars more massive than the sun. 



line of Fig.[3]for the hottest star shows the AQ values slightly 
exceeding the largest d ifferential rotation observed to date 
jjeffers fc Donatill2008l ). 

We used statistics of the results of Fig. [2] and [3] and also 
the results of computations for stellar masses intermediate 
between the lines of these Figures to define the coefficients 
of the regression (|19|l . Only the results for rotation periods 
from slightly less than 1 day to P™ t ax of Eq. (f2"U)l were in- 
cluded. Total statistics consists of about 1600 entries from 
differential rotation computations. The neural network ()19p 
with N = 6 gives quite a satisfactory approximation with a 
typical error of about 3%. Coefficients of the regression are 
given in Table 1. 

Given the observed rotation period and temperature of 
a star, a theoretical prediction for the surface differential 
rotation can easily be estimated by using coefficients of this 
Table. It should be noted that the estimation applies to sin- 
gle main-sequence dwarfs only. It cannot be used for giants, 
close binaries, or pre-main sequence stars. 



Table 1. Coefficients of the neural network approximation H19H 
for the surface differential rotation in terms of the effective tem- 
perature and rotation period. 
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Figure 4. Angular velocity isolines in the convection zone and 
the angular velocity profile on the surface of IMq stars rotating 
with periods of 30 days (a) and 3 days (b). 



For a given mass, the amount of the surface differential 
rotation of Fig. [2] and [3] varies by about 30% as the rotation 
period changes. This mild variation is, however, worth dis- 
cussing. The surface differential rotation increases initially 
up to a maximum but then decreases with P ro t- Maxima on 
all lines are positioned at the values of the Coriolis number 
(0 between 10 and 20 (the Coriolis number is estimated for 
the middle of convection zone). Stars of smaller mass host 
slower convection with larger r and attain their largest dif- 
ferential rotation at smaller rotation rates. The increase of 
differential rotation with angular velocity for not too rapid 
rotati on was also found in 3D simulations of iBrown et al.l 
(2008). The 3D simulations predict, however, much larger 
differential rotation compared to our mean-field model. 

The character of differential rotation changes between 
cases of slow (fl* < 10) and rapid (S7* > 20) rotation. Fig- 
ure [4] shows the computed internal rotation of IMq stars 
rotating with periods of 30 and 3 days. The slow rotator 
has a smooth distribution of angular velocity on the sur- 
face and inside its convection zone. Isorotational surfaces 
are far from cylinders. The rapid rotator has thin bound- 
ary layers near the top and bottom of its convection zone. 
The layers result from violation of the Taylor-Proudman bal- 
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Figure 5. Meridional flow stream lines and the flow velocity as 
function of radius at the latitude of 45° for IMq stars rotating 
with periods of 30 days (a) and 3 days (b). Negative velocity 
means poleward flow. 



ance near the stress-free boundaries. The boundary layer s 
wer e found in mean-field c o mput ations of IPurnevI (jl989T ) 



and iKitchatinov fc Riidigerl dl999l ) as well as in 3D simu 

it 



lations of iBrun fc Toomrel (|2002l ) and iBrown et al.l |200St ). 
The profile of angular velocity on the surface of the rapid 
rotator has a 'peculiarity' around the latitude where the an- 
gular velocity isoline tangential to the inner boundary at the 
equator arrives at the surface. Isorotational surfaces for the 
rapid rotator are much closer to cylinders. 

Isorotational surfaces are cylinder-shaped near the 
equator and disk-shaped near the axis of rotation in all our 
computations. Isorotational surfaces are normal to the equa- 
torial plane due to the symmetry of angular velocity dis- 
tribution about the equator. As the distribution should be 
regular at the poles, isorotational surfaces are also normal 
to the rotation axis. Therefore, it is a general, though ele- 
mentary rule, that isorotational surfaces are cylinder-shaped 
near the equator and disk-shaped near the poles (for an al- 
ternative opinion see iBalbus! [20091 ). 

Violation of the Taylor-Proudman balance in the 
boundary layers results in relatively large values of the 
sources of meridional flow written on the right side of (|10p . 
Accordingly, the flow attains its largest velocities close to 
the boundaries. Figure [5] shows the meridional flow struc- 
ture for the same computations as Fig. 2] The boundary lay- 
ers in slow rotators are relatively thick (KOI I). The flow is 
smoothly distributed over the entire thickness of the convec- 
tion zone in this case. In a rapid rotator, however, the flow 
is confined in thin lay ers near the boundaries. This a grees 
with 3D simulations of I Brown et al. | (|2008|j2010l . l201lh who 
found a decrease of the meridional flow energy with rotation 
rate. 

The near-boundary flows are significant. The surface 
flow is important because it is potentially observable. The 
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Figure 6. Dependence of meridional velocities on the top (full 
lines) and bottom of convection zone (dashed) on rotation period 
for stars of different mass. The velocities are taken at the latitude 
of 45° . Negative values mean poleward flow. 



flow near the bottom of the convecti on zone is increa singly 
recognized as important for dynamos (|Choudhuril201ll ) . Fig- 
ure [6] shows how the near-boundary flows vary with stel- 
lar mass and rotation rate in our computations. The bot- 
tom flow is smaller but not much smaller compared to the 
surface. Note, however, that the flow does not penetrate 
deep beneath the convection zone. The meridional veloc- 
ity decreases rapidly wi th depth below the convection zone 
( Oilman fc Mieschll2004l ). Similar to differential rotation, the 
flow amplitude increases with stellar mass and it increases 
steadily with rotation rate. The flow in the bulk of the con- 
vection zone, however, has the opposite tendency to decrease 
with rotation rate (Fig. [5| . 

The meridional flow is produced, in particular, by the 
effect of thermal wind due to the temperature dependence on 
latitude. Figure[7]shows the depth profiles of the 'differential 
temperature' 



ST(r) = - (S(r,e) \e=o -S(r,6) | 9=} ) 



(21) 



for the same computations as Figs. [4] and [5] The differen- 
tial temperature in mean-field models is produced partly by 
the dependence of thermal eddy conductivity on latitude but 
mainly by anisotropy of the thermal conductivity tensor ()13|) 
llRiidiger et al.ll2005l ). The anisotropy is induced by rotation 
and it increases with rotation rate. Accordingly, the differ- 
ential temperature in faster rotating star of Fig. [7] is larger. 
There were many attempts to observe the small tempera- 
ture difference betwee n the equator and poles of the sun. 
Recent observations of iRast et alj l|200ct ) suggest that the 
solar poles are warmer than the equator by about 2.5 K. 

The anisotropy of thermal eddy conductivity is very im- 
portant for differential rotation formation. It was not pos- 
sible to reproduce the differential rotation of the sun ne- 
glecting this anisotropy (|Riidiger et al.1 12005). Reducing the 
anisotropy, e.g., twice results in a change in the sense of 
differential rotation to an anti-solar state (the equator ro- 
tating slower than the poles) in slowly rotating stars. Only 
in this way were we able to reproduce a similar change in 
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Figure 7. Entropy isolines and depth profiles of the pole-equator 
temperature difference ( 12 1 D for IMq stars rotating with periods 
of 30 days (a) and 3 days (b). 
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Figure 8. Differential rotation of 0.5Mq (a) and 1.2Mq (b) stars, 
both with P ro t = 10 days. 



th e sense of different ial ro tation found in th e 3D s imulations 
of iMatt et all i|201ll ) and lBessolaz fc Brunl l|201ll) . 

Whether the differential rotation of a star belongs to the 
slow or rapid rotation regime is not solely controlled by the 
rotation period. Figure [8] compares the rotational states of 
O.5M0 and 1.2Mq stars with the same P ro t = 10 days. The 
more massive star with Q* — 3.7 is a typical slow rotator 
while the smaller star with f2* = 49 is not. 



4 SUMMARY 

Our main purpose was to predict the surface differential ro- 
tation in terms of measurable stellar parameters. This pur- 
pose is achieved by formulating the analytical approximation 
(|19fl for the differential rotation of main-sequence dwarfs as 
a function of their surface temperature and rotation period 
(Table 1 gives the coefficients of the approximation) . 

The prediction is based on the mean-field model of dif- 
ferential rotation (KOll), which relies on the quasi-linear 
theory of turbulent transpor t coefficients for rotating fluids 
l|Riidiger fc Hollerbacdl20olh . The results of extensive com- 
putations for stars of d ifferent mass and ag es were processed 
using Gyrochronology |Barnesll2003l 120071 ) and temperature 
scaling for differential rotation. The scaling means that the 
dependence of differential rotation on chemical composition 
and mass can be combined into a common dependence on 
temperature. The temperature scaling and Gyrochronology 
help to reduce the number of stellar parameters on which 
the differential rotation depends to just two - rotation pe- 
riod and temperature. 

Our computations suggest that the main trend in the 
dependence of differential rotation on stellar parameters is 
its increase with effective temperature. The increase steep- 
ens for F-stars compared to cooler stars (Fig. [TJ so that the 
maximum surface differential rotation, AS! ~ 0.6radday~ , 
is achieved for the hottest F-stars we considered. There is 
also a weaker dependence on rotation rate, which is not 
monotonous. The surface differential rotation of young stars 
increases initially by about 30% as the star ages and rota- 
tion period increases but then it changes to a decrease with 
Prot- The maximum differential rotation is achieved at the 
values of the Coriolis number (0 between 10 and 20. This 
range of the Coriolis number separates the regimes of fast 
and slow rotation. Angular velocity distribution in slow ro- 
tators is smooth. Rapid rotators have thin boundary layers 
near the top and bottom of their convection zones where 
angular velocity and meridional flow vary sharply. 

One-cell meridional flow with poleward flow on the sur- 
face and a return flow at the bottom of the convection zone 
was found in all our computations. The flow attains its max- 
imum velocity on the top boundary. The amplitude of the 
flow increases with the effective temperature, but not as 
much as differential rotation does, and remains of an or- 
der of several tens of meters per second. The flow ampli- 
tude increases smoothly with rotation rate. The structure 
of the flow, however, changes considerably between cases of 
slow and fast rotation. The flow is distributed smoothly over 
depth in the convection zone in slow rotators. In rapid ro- 
tators, however, the flow is concentrated in the boundary 
layers near top and bottom with very weak meridional cir- 
culation in the bulk of the convection zone. 
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